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Оценка параметров распределения 
дефектов ИС с помощью 
помехоустойчивой кластеризации 


Проведена оценка параметров распределения дефектов ИС с помощью помехоустойчивой кластеризации. 
Предложены метод и процедура оценки размеров блока среднекластерного распределения дефектов 
ИС с использованием мультистартового субградиентного итеративного метода кластеризации. Оценена 
помехоустойчивость метода кластеризации и его погрешность. 


Развитие микроэлектроники требует обеспечить повышение технологического 
уровня и экономической эффективности полупроводникового производства. При этом 
возрастает актуальность формирования адекватных статистических распределений, 
позволяющих на стадии анализа технологического процесса производства прогнози- 
ровать распределение дефектов, которые вызывают ошибки функционирования схемы, 
и выход годных ИС. Поэтому при производстве устойчивых к наличию дефектов ИС 
более точное определение выхода годных определяет количество добавляемых к схеме 
избыточных элементов. Такой подход характерен для схем, состоящих из модулей, 
таких как элементы памяти и цифровой логики. Для прогноза выхода годных в на- 
стоящее время применяются различные варианты компаунд-распределений на основе 
обобщенного распределения Пуассона [1], [2]. Компаунд-распределение позволяет 
учитывать неоднородность распределения дефектов, приводящих к потере годных 
кристаллов в процессе производства ИС. Эти дефекты не распределены внутри пла- 
стины равномерно, а образуют кластеры [1]. В работах [1], [2] обоснован иерархический 
подход к оценке выхода годных ИС. Эта иерархия включает в себя крупнокластер- 
ное, среднекластерное и мелкокластерное распределения. Выбор вида распределения 
обусловлен соотношением размеров кластера, который характеризует неоднородность 
распределения дефектов в анализируемой области, и самой анализируемой области 
(модуль на полупроводниковой пластине, пластина, набор пластин (партия)). Понятие 
кластера в работах не определено [2]. В работе [1] введено понятие «блока», предпо- 
лагается, что дефекты коррелированы между собой только внутри этих блоков, сами же 
блоки статистически не зависимы друг от друга. Выбор размеров блоков проводится 
с учетом типа ИС, особенностей производственных линий и процессов обработки 
пластин [1]. 

Ошибки оператора, сбои оборудования при измерении обусловливают зашум- 
ленность данных о дефектах при малом объеме их выборок. Методы, основанные на 
сравнении распределений этих данных, и распределения на их основе чувствительны 
к таким отклонениям (шуму) данных [1]. Влияние этих недостатков может быть умень- 
шено при применении помехоустойчивых процедур кластеризации. В связи с этим 
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представляется актуальным разработать метод, позволяющий формализовать выбор 
размера блока среднекластерного распределения оценки выхода годных, применяя 
помехоустойчивые методы кластеризации, принятые при распознавании образов [3]. 

При кластеризации данные разделяются на кластеры по признаку компактно- 
сти, так, чтобы был оптимизирован функционал качества. Метод оптимизации выбирают 
в зависимости от особенностей формирования и свойств этого функционала качества, 
который может быть явно не известен, может обладать поверхностью многоэкстре- 
мальной, зашумленной, поскольку анализ производится по малым выборкам. Сущест- 
вующие методы кластеризации в этих условиях работают плохо. Для снижения влияния 
указанных выше недостатков в работе предлагается применить разработанный авто- 
рами помехоустойчивый субградиентный итеративный метод кластеризации в прост- 
ранстве вейвлет-преобразования [3]. 

Целью данной статьи является разработка метода, позволяющего формализо- 
вать расчет размеров блока среднекластерного распределения оценки выхода годных 
с применением помехоустойчивой кластеризации. Для достижения этой цели реше- 
ны задачи анализа методов оценки выхода годных в микроэлектронике; разработки 
процедуры оценки выхода годных и размеров блока среднекластерного распределения с 
использованием субградиентного итеративного метода кластеризации; проведения экс- 
периментальных исследований при оценке помехоустойчивости метода кластеризации. 

При кластеризации используют иерархические и итеративные методы. Обе группы 
методов чувствительны к шуму данных [3], поэтому целесообразно разрабатывать 
методы, снижающие чувствительность к шуму при кластеризации. 


Анализ методов оценки выхода годных ИС 


При оценке выхода годных ИС используют компаунд-распределения дефект- 
ности [1], [2]. Дефекты, приводящие к потере годных при производстве ИС, не 
распределены внутри пластины равномерно, а образуют кластеры. Вероятность иметь 
в выборке К дефектных изделий определяется обобщенным отрицательным биномиаль- 
ным распределением (компаунд-распределением) [2]. 


СТА ых _ Га)“ 
Ро - и° 2 = тазьее 


гдеа> 0,6 >0; Г(а) — гамма-функция; а — параметр формы; Ъ — параметр масштаба. 
Выход годных (вероятность получить К = 0 дефектных модулей) в рамках ком- 
паунд-распределения Пуассона для крупнокластерного распределения, когда вся 
микросхема располагается внутри одного блока, определяется 
1 
> 
(1 + МА 


ш 


УЕР(К=0) = 


где ^„ - среднее число дефектов на один модуль; @а„ — параметр кластеризации; 
В- размер блока; М- размер микросхемы [1], [2]. 
Для среднекластерного распределения, когда 1 < В < М, выход годных модулей 


У=Р(К = 0) = 


450 «Искусственный интеллект» 372009 


Оценка параметров распределения дефектов ИС... 8 щ 


Для мелкокластерного распределения выход годных модулей 


1 


(+ =) № 
(7 


ш 


У=Р(К = 0) = 


Среднекластерное распределение включает в себя крупнокластерное и мелко- 
кластерное распределения, как предельные случаи [1]. Размер «блока» для мелко- 
кластерного распределения равен одному модулю, для крупнокластерного распреде- 
ления — равен всей площади пластины в модулях. Среднекластерное распределение 
применяют для микросхем большой площади (порядка квадратного дюйма), когда 
размер микросхемы больше размера кластера. В [1] введено понятие «блока» как еди- 
ницы измерения размера кластера. Блоки статистически не зависимы друг от друга, 
внугри них дефекты распределены равномерно и коррелированы между собой, число 
дефектов в блоке имеет отрицательное биномиальное распределение [1]. Коэффици- 
ент кластеризации © постоянен внутри блока, но © возрастает, когда рассматривае- 
мая область включает в себя несколько блоков [1], [2]. Для определения размеров 
блока разработан ряд методов. При определении размера блока согласно первому ме- 
тоду сначала оцениваются коэффициенты кластеризации © для каждого потенциального 
блока размером 1х 1, 1х2, ...1,7] модулей. Из этих значений о(Т,Т) формируется мат- 


рица. Далее проводится поиск наибольших (1,7) с (ГУ) меньшим, чем @(1). Это 
количество модулей принимается как размер блока [1]. Во втором методе оценки 
размеров блока размер блока определяется в два этапа (сначала количество модулей 
в длину В, , потом — в ширину В,). Пластина делится на блоки размером (1,7) моду- 


лей. Чтобы определить В, и В, ‚ формируются две матрицы с результатами проверки 


по критерию х? статистической независимости между соседствующими вертикаль- 


но (горизонтально) блоками. Индексы Ги У, которые определяют выбор В; и В,, 


определяют по первой строке (столбцу) матриц, в которых значения х” ниже, чем в 
других. Однако в работе [1] те же авторы констатируют, что, даже применяя оба метода 
в сочетании, определить размер блока в случае малых выборок (малого количества 
исследуемых пластин) сложно. В работе [1] при определении параметров среднеклас- 
терной модели пластину при производстве ИС делят на модули, формируя из них 
«блоки» размером 1х 1, 1х2, ..., 2х З и так далее модулей, определяют для каждого ва- 
рианта параметры распределения с, /. ‚ размер блока. Выбирают тот вариант, который 


по критерию Колмогорова дает наименьшее расхождение с результатами эксперимента. 

Размер блока, который используют для определения коэффициента кластеризации, 
итеративно определяется заново для каждого вида производства ИС [1]. Это опреде- 
ление в пределах пластины предполагает перебор различных сочетаний количества и 
расположения модулей, входящих в блок. Все варианты распределений выхода год- 
ных ИС чувствительны к шумовым измерениям при контроле по малым выборкам [1]. 
Поэтому представляется актуальным предложить метод, позволяющий формализовать 
и автоматизировать процесс определения размеров блока для прогноза выхода год- 
ных ИС с помощью среднекластерного распределения на базе помехоустойчивого 
мультистартового субградиентного итеративного метода кластеризации в пространст- 
ве вейвлет-преобразования [3]. 
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Мультистартовый метод кластеризации 


Задача кластеризации состоит в разбиении множества образов объектов на группы 
(кластеры) с учетом их сходства, которое определяют через расстояние. Расстояние в 
выбранном подходе к кластеризации рассчитывается от объектов к центру кластера. 
Координаты центров кластеров находятся одновременно с разбиением данных на 
кластеры. При отсутствии априорной информации о форме исследуемых кластеров и 
с учетом предполагаемой зашумленности выборки принимается функционал из группы 
родственных функционалу минимума дисперсии [3]. 

Для решения задач оптимизации при таких условиях авторами разработан по- 
мехоустойчивый субградиентный итеративный метод оптимизации в пространстве 
вейвлет-преобразования [4]. Для повышения помехоустойчивости кластеризации на 
основе субградиентного итеративного метода оптимизации авторами разработан мульти- 
стартовый субградиентный итеративный метод кластеризации в пространстве вейвлет- 
преобразования [3]. 


При итеративном подходе к кластеризации определяют вектор с =с,„, который 


оре > 
доставлял бы экстремальное значение О(х, ©) — функционалу вектора переменных 
с= (с, ,...ск), зависящему от вектора случайных последовательностей х=(х/,...хм) [3]. 


По показам образов х Е Х определяют центры множеств Х,, и их границы. При этом 
м 

О(х, с1,.... См) = Ук (Х, С1,..., См )Ех (х, с1,....См) -— реализация функционала ка- 
К=1 


чества; Е, (х, с|,...,См)- функция расстояния элементов х множества Х от центров 
с, кластеров Ху; &, (:) - характеристические функции 
1, когда хе Х,, 


ЕО Сами ба = 
кс м) 0, когдахя Х,. 


Поисковый итеративный алгоритм для определения центров двух кластеров с" и с> 
с [и] = с [и -П-у[иМ „Об, с [и — Ц, с,[п 1) 
с, [п] = с[и-П-у[и № „„„ Оса, с[и-— Ц, [и 1), 
где У. ОС], < [п-П, с, ш-1), У „О, < [в -П, с, [п -1]) — оценки градиента 
реализации для первого и второго кластеров соответственно; у„[п] — шаг; п — номер 
итерации; К — номер кластера. 

Но, если в условиях помех оценка градиента проводится разностным методом, 
поисковый регулярный итеративный метод при кластеризации даст низкую помехо- 
устойчивость. 

Метод кластеризации в пространстве вейвлет-преобразования заключается в 


следующем. После инициализации параметров метода, при кластеризации в каждый 
момент времени для каждого из 1 элементов взвешенной суммы с вейвлет-преобра- 


зованием определяют значение характеристических функций 51 (х,с|,с.),1=1,2 ‚ вхо- 
дящих в оценку субградиента. Для этого по методике [3] пары значений 

с [в-П, с [в-П; си -Пюва[п], с›[п-П; с [п -П, с›[п- | +1еа[п] (1=ь М) 
при данном х[п] подставляют в Х(х, с, со) = |х[п] - с1 | - [п] - с2 Г. Здесь М — длина 


носителя вейвлет-функции; а[п] — скаляр. Функция Ё(х, с,,с,) равна нулю на границе 
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и имеет различные знаки в различных областях. Поэтому, если Ё(х,с,,с>) отрица- 
тельна, = =1,=, =0, если положительна, &, =0,&, =1 [3]. 

В качестве базового использован градиентный метод поиска минимума [3], [4]. 
Исходные данные для его работы: начальные значения координаты минимума и шага у = 1, 
коэффициент, обусловливающий изменение шага ‘у вблизи минимума В = 0,5, точность 
определения оценки градиента =, количество итераций ]. Процедура вычисления ми- 
нимума при кластеризации включает [3]: 

— вычисление оценки градиента; 
— если значение оценки градиента меньше заданного значения точности 5 — останов; 
— вычисление величины шага: задается начальное значение величины шага у=1; 


вычисляется вспомогательное значение приращения функции Д, если приращение 
функции Д меньше нуля — у[п] =у и переход к следующему этапу, иначе [п] = Ву 
и переход к предыдущему этапу; 
— расчет координаты минимума на п-ой итерации; 
— п=0+1 и переход к начальному этапу вычисления минимума при кластеризации. 

При вычислении оценки градиента на каждой итерации на первом этапе вычисляется 
взвешенная сумма значений минимизируемого функционала О(х[п], с [п-1, с [п-1]) 
с вейвлет-функцией Хаара. Это позволяет переместить поиск в район экстремума с 
погрешностью, определяемой асимметрией этого функционала. На втором этапе оцен- 
ки градиента при кластеризации вычисляется взвешенная сумма минимизируемого 
функционала О(х[1], с [п-П, с, [п-1]) с гиперболической функцией (1) = Е при 

ох 
начальном масштабе са = 0,5. 
НУТ(<]) = Ох], с [п — Ц, [п -1)*Ф 0, 

где * — операция взвешенного суммирования. 

Далее определяют приближение к значению координаты центра кластера по 
итеративному алгоритму в пространстве вейвлет-преобразования по схеме 

с +П=с а] + у Н\Т (<), 

где Н\Т(с[п]) — взвешенная сумма с вейвлет-функцией в точке с[п]; У[п] — шаг. 

Если найденная на этом этапе координата оптимума отличается от найденной 
на предыдущем этапе не более, чем на 6 (заданную точность поиска координаты оп- 
тимума), поиск заканчивается. Для оценки субградиента использовано гиперболи- 
ческое вейвлет-преобразование, полученное по следующей схеме. На каждом уровне 
поиска координаты оптимума значение масштаба © увеличивается в соответствии с 
© = {0,5; 1; 2; 3; 4; 5}. Если условие окончания поиска координаты оптимума при зна- 


чении величины 0 =5 не достигается, оценка субградиента производится разностным 
методом. После этого поиск заканчивается. При поиске координаты оптимума осущест- 
вляется переход от поиска с помощью вейвлета Хаара, с высокой помехоустойчи- 
востью, вплоть (с ростом %) до поиска с помощью дифференциатора, способного дать 
максимальную точность. Далее проверяют условие точности определения центра клас- 
тера, если он достигнут — останов. 

Оценка помехоустойчивости метода кластеризации проводилась с использованием 


функции #(х)=х” при значениях ее аргумента х =1,...80. Помеха была распреде- 


лена по нормальному закону с нулевым средним и среднеквадратическим отклонением, 
изменявшимся в диапазоне от 0 до 5477. При отношении сигнал/Апум по амплитуде 1,17 
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(помеха распределена по нормальному закону с нулевым средним и среднеквадрати- 
ческим отклонением 5477, максимальное значение тестируемой функции Г(х) = 6400) 


относительная погрешность определения минимума составила 8,32% [3], [4]. 


Оценка параметров среднекластерного распределения 


Для оценки параметров среднекластерного распределения предлагается следую- 
щая процедура. 

1. Разделение полупроводниковой пластины на модули равной площади [1]. 

2. Определение Д„,@„ — количество дефектов на модуле, коэффициента клас- 


т 


теризации для модуля соответственно проводится по методике [2]. Тогда 


где \ — оценка /,_; $ — количество пластин; \! — количество модулей на пластине; 
ш ш 2 > 


Х; -— число дефектов на модуле 1 (1=1....,5: \); Х - среднее число дефектов на 
модуле. 
Оценка ©, проводится с учетом того, что число дефектов Х имеет отрицатель- 


ное биномиальное распределение с параметрами Д„, &„. Оценка дисперсии будет 


т? 


УС® = 1 Ух ве 
5. М 
Оценка коэффициента кластеризации для модуля 
Х2 
— УСО-Х. 
Коэффициент кластеризации блока определяется как @, = @ и [1]. 


бы 


3. Определение количества, центра, диаметра кластеров с заданной точностью 
9 с помощью предложенного метода кластеризации. 

4. Определение максимального целого числа модулей, входящих в кластер, для 
каждого варианта кластеров, определенных в пункте 3. 

5. Формирование вариационного ряда площадей кластеров, выраженных в мо- 
дулях. Эти значения (1,7) принимаются как значения размеров блоков для первой 
итерации. 

6. Определение размеров блока. 

6.1. По методике [1] последовательно 
— формируется матрица оценок коэффициентов кластеризации © для потенциальных 
размеров блока от 1Ж1 и до /х./ модулей; отыскивается наибольшее (1,7), для кото- 


рого @(Т.Т) меньше, чем 0(1) , это сочетание количества модулей и принимается как 
размер блока; 

— применяется второй метод оценки размеров блока; размер блока этим методом 
определяется в два этапа (количество модулей в длину и в ширину соответственно 


В, иВ, ). Исследуемая пластина делится на блоки размером (1,7) модулей. Чтобы оп- 
ределить В| и В, ‚ формируются две матрицы, содержащие результат проверки по кри- 
терию х2 статистической независимости между соседними блоками. Индексы [и Т, 
которые определяют выбор В| и В, ‚ определяются по первой строке и столбцу со- 


ответствующих матриц, в которых значения х2 значительно ниже, чем в других. 
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6.2. Последовательно, начиная с минимального (1,7), сравниваются значения (1,7) 


из вариационного ряда, полученного в пунктах 5 и 6.1. При совпадении результатов с 
заданной точностью ф, принимается минимальный размер блока, для последующего 


сравнения с результатами экспериментальной оценки выхода годных. 
Если оценки размеров блока в пунктах 5 и 6.1 отличаются более чем на ф, ре- 


комендуется повысить точность кластеризации (уменьшить 0) и повторить расчеты 
по пунктам 3 -— 6. Если различаются все три оценки размеров блока по пунктам 5 и 6.1, 
выбирают тот вариант размера блока, который при сравнении с результатами экспе- 
риментальной оценки выхода годных дает меньшую ошибку. 

7. Определяют количество дефектов блока № = В - А». Здесь В = Ву: В› — размер блока. 

В работе решены задачи анализа методов оценки выхода годных в микроэлек- 
тронике, разработки метода и процедуры оценки размеров блока среднекластерного 
распределения дефектов ИС с использованием мультистартового субградиентного итер- 
ативного метода кластеризации; проведены экспериментальные исследования по оценке 
помехоустойчивости метода кластеризации и снижения его погрешности. Относитель- 
ная погрешность определения минимума при кластеризации для тестовой функции при 
отношении сигнал/шум по амплитуде 1,17 составила 8,32%. Эти результаты позво- 
ляют рекомендовать разработанный метод определения размера блока среднекластер- 
ного распределения в процессе статистического анализа годных ИС при высоком 
уровне помех и при малых объемах выборок. 
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Ощнка параметрив розпод1лу дефектв ГС з допомогою завадост ко! кластеризаци 

Проведено ощнку параметр розподму дефект!в [С з допомогою завадостйко! кластеризаци. Запропонован! 
метод 1 процедура ощнки розм!ртв блоку середньокластерного розподлу дефектив ГС з використанням 
мультистартового субградентного 1теративного методу кластеризаци. Проведена ощнка завадостйкости 
методу кластеризаци та його похибки. 


С. Уи. 5йсйегтракКоъа, Г.М. Кгуоу, Г.С. Аракитоу 

ТС ОеГесе Оба ийоп Рагашеег$ Еуашайоп ууйВ Мое За Шу Сшяегто 

ТБе 1С деесё ЧзитБийоп рагате{ег$ еуашайоп УИ по15е заб сГазегие уаз сатутеа оп. ТВе ТС 
тет сГауегие 4еесё 45а Ноп?”5 ЫоскК $12е еуааНоп тефо4 апа 15 паретешщайоп ргоседиге \аз 
ргорозеа. шп фай ргоседиге то] $агипе заб огаЧеп{ Иегайуе стазегте те о4 15$ изе4. ТВе ехрегипета1 
шуезйраноп Юг 1$ сГазегте тефо4 по15е плтипйу шсгеазште ап4 о\уп еггог десгеазте езйтайоп ууаз 
сагле4 ой. 
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